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l~:> ' Abstract. Good understanding of the cross sections for (anti)neutrino scattering off nuclear targets in the few-GeV energy 

^Nj , region is a prerequisite for the correct inteipretation of results of ongoing and planned oscillation experiments. To clarify 

5, , a possible source of disagreement between recent measurements of the cross sections on carbon, we analyze the available data 

O ■ within an approach based on the realistic spectral function of carbon, treating neutral-current elastic (NCE) and charged- 

^rfH ' current quasielastic (CCQE) processes on equal footing. We show that the axial mass from the shape analysis of the 

^^ I MiniBooNE data is in good agreement with the results reported by the BNL E734 and NOMAD Collaborations. However, the 

[~^ . combined analysis of the NCE and CCQE data does not seem to support the contribution of multinucleon final states being 
large enough to explain the normalization of the MiniBooNE-reported cross sections. 

Keywords: neutrino-nucleus scattering, charged-current neutrino interactions, neutral-current neutrino interactions 
'^ '. PACS: 13.15.+g, 25.30.Pt 

O^- INTRODUCTION 

The correct interpretation of the outcome of oscillation experiments requires a precise knowledge of the (anti)neu- 
trino cross sections. This is the case even in those experiments in which the event yields in near and far detectors are 
compared, because, in general, flux and backgrounds do not scale in a simple manner [1]. 

The description of nuclear effects in neutrino scattering is now generally regarded as one of the main sources of sys- 
tematic uncertainties in oscillation experiments. In particular, profound understanding of charged-current quasielastic 
[■-«. ' (CCQE) interaction with nucleons bound in the nucleus plays a crucial role. 

t;J" . The carbon nucleus is of special importance, because the targets applied in neutrino detectors often involve carbon 

._b ' compounds, such as mineral oil, polystyrene, or organic scintillators. The total 'gC(v^,/i^) cross section as a function 

f^ , of energy has recently been measured with the NOMAD [2] and MiniBooNE [1] expeiiments. 

C^ ' In NOMAD, studying neutrinos of energy down to 3 GeV, CCQE events with and without knocked-out proton 

detected have been analyzed separately and, after adjusting the description of final-state interactions, shown to yield 
consistent results. Thanks to the 45-GeV average energy of the deep-inelastic-scattering events, the normalization has 
been determined from the well-known total inclusive charged-current (CC) cross section and from the purely leptonic 
process of inverse muon decay. The observed reduction of the cross section due to nuclear effects is ~4%. 
^ I In MiniBooNE, all events without pions detected have been classified as CCQE and the total cross section is 

■ ■ ■ extracted for energy up to 2 GeV. The high statistics of CCQE events (^ 1 x those in NOMAD) allowed for obtaining, 

for the first time, the double differential cross section. Surprisingly, the repotted total cross section for carbon is higher 
than that for free neutrons. The size of the effect is ^5% at a few hundred MeV, increasing to ^15% at neutrino energy 
higher than -850 MeV. 

As a 20% uncertainty of the cross section would have an impoitant impact on determination of oscillation parame- 
ters [3], the difference between the NOMAD and MiniBooNE results requires very careful theoretical analysis. 

In this paper, we cover both NCE and CCQE processes, treating them in an identical manner and comparing the 
obtained cross sections to those measured by the Brookhaven National Laboratory Experiment 734 (BNL E734) [5], 
MiniBooNE [1,6], and NOMAD [2]. Our analysis is presented in detail in Ref. [7]. 

DESCRIPTION OF THE APPROACH 

We assume that the process of neutrino-nucleus interaction involves a single nucleon and the remaining (A — 1) 
nucleons act as a spectator system. In this regime, called the impulse approximation (lA), the neutrino-nucleus cross 
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FIGURE 1. Differential NCE cross sections a^'^/dQ^ in (a) BNL E734 and (b) MiniBooNE. The SF calculations are compared 
with the data from Refs. [5] and [6], respectively. The error bars show the statistical and systematic uncertainties added in 
quadrature. In panel (b), the theoretical result multiplied by 1.2 is also shown for comparison. 

section is obtained convoluting the elementary neutrino-nucleon cross section with the hole and particle spectral 
functions (SFs) of the nucleus. The hole spectral function Phoie(Pi£^) is the probability distribution of removing 
a nucleon of momentum p and leaving the residual nucleus with excitation E in its rest frame. The particle spectral 
function fpart(p', ■'^') describes the propagation of a nucleon of momentum p' and kinetic energy ,'^' . 

Realistic hole spectral functions for various nuclei have been obtained by the authors of Ref. [8] in the local density 
approximation (LDA), consistently combining the shell structure of the nucleus with the correlation contribution 
obtained from theoretical calculations for uniform nuclear matter at different densities [8, 9]. The carbon SF of Ref. [8], 
employed in this paper, has been extensively used in the analysis of electron scattering data in various kinematical 
setups. We use the particle SF obtained in the LDA scheme by convoluting the momentum distribution of nuclear 
matter at different densities with the density profile of the carbon nucleus [4]. 

To account for multinucleon processes, in scattering on nucleons bound in the carbon nucleus, the effective axial 
mass Ma = 1 .23 GeV is applied. While this purely phenomenological method cannot be expected to be accurate in 
any kinematical regime, its validity is quantitatively supported by the results of Nieves et al. [10] for the double 
differential cross section in a broad kinematical range of MiniBooNE. Moreover, in the context of our work, the 
kinematical setup of NOMAD does not differ significantly from that of MiniBooNE. It is a simple consequence of the 
fact that (quasi)elastic processes constrain the high-|q| contribution to appear solely at high Q^, making it negligibly 
small due to the nucleon form factors. Hence, the phenomenological approach to multinucleon processes seems to be 
applicable also in comparisons to the NOMAD results. 

In this paper, we consider the total cross sections and the flux-averaged differential cross sections da jdQ^ for 
neutrino and antineutrino (quasi)elastic scattering. As FSI may cause only a redistribution and a shift of the strength, 
they do not affect the total cross section. The differential cross sections da jdQ^ are expected to be modified by FSI 
at Q^ < 0. 15 GeV^ [11], because only at this kinematics the real part of the optical potential significantly changes the 
typical energy of the knocked-out nucleon. However, in the low-Q^ regime, the validity of the impulse approximation, 
underlying our calculations, becomes questionable [4]. Therefore, in this article, FSI are not taken into account. 



COMPARISON TO THE NCE AND CCQE DATA 

From the neutrino experiments which obtained results on NCE scattering off nuclei, the highest event statistics to date 
have been collected by the BNL E734 and MiniBooNE Collaborations. 

In the BNL E734 experiment [5], studying vp and Vp NCE interactions, the target was composed in 79% of protons 
bound in carbon and aluminium and in 21% of free protons. The mean value of neutrino (antineutrino) energy was 
1.3 GeV (1.2 GeV). Determination of the v and v fluxes involved fitting to the CCQE event sample. Based on 
1686 (1821) candidate vp (vp) events surviving the cuts, the flux-averaged differential cross sections da^'" /dQ^ 
for scattering off the BNL E734 target were extracted [5]. 

Figure 1(a) shows that the SF calculations provide a fairly good description of the BNL E734 data [5] in the whole 
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FIGURE 2. Total CCQE cross sections for muon (anti)neutrino scattering off carbon. The SF results are compared with the data 
extracted from the (a) MiniBooNE and (b) NOMAD experiments. The error bars show the total uncertainties. In panel (a), the 
calculation multiplied by 1.2 is also shown for comparison. 







// V,,, MiniBooNE ^^^ - 
// SFxl.2 




range of Q^. Note that the agreement seems to be better with the higher-statistics antineutrino data. This is also the 
case for neutiinos in the region of the lowest uncertainty, 0.5 < Q^ < 0.8 GeV^. 

The MiniBooNE experiment, using the Cherenkov detector filled with mineral oil (CH2), has been sensitive to 
both vp and vn NCE scattering [6]. In neutrino mode, its beam has been composed almost exclusively of muon 
neutiinos, with an average energy of 788 MeV. The neutiino flux at the detector has been determined by a Monte 
Carlo simulation [12]. The accuracy of the shape determination has been proven by the comparison of the observed and 
predicted v^ energy distribution in the CCQE event sample. However, the normalization of the measured distribution 
of CCQE events is reported to be higher by a factor of 1.21 ±0.24 than the calculated one [12]. The MiniBooNE 
Collaboration recorded in neutrino mode 9453 1 candidate events surviving the NCE cuts, which allowed extraction of 
the differential cross section with unprecedented precision. 

Figure 1 (b) shows that our calculations reproduce the shape of the MiniBooNE NCE v cross section very accurately, 
both in the region of lower ^ and in the tail. However, to match the absolute scale of the MiniBooNE data, it was 
necessary to multiply the SF results by a factor of 1.2, consistent with the normalization discrepancy of 1.21 ±0.24 
observed by the MiniBooNE Collaboration in the CCQE analysis [13]. 

To clarify a possible source of the normalization discrepancy, it is useful to check how our approach describes 
the CCQE data. From the collected in neutrino mode 146070 events passing the cut, the MiniBooNE Collaboration 
has extracted, among others, the flux-unfolded total CCQE cross section [1]. As presented in Fig. 2(a), the energy 
dependence of our result is in good agreement with the experimental points, however, also in this case the factor of 1 .2 
is required to match the normalization. 

From our analysis, it consistently emerges that, while correctly describing the shape of the MiniBooNE-reported 
cross sections, our calculations fail to reproduce their absolute scale, remaining inadequate by 20%. 

In this context, it is interesting to compare our results to the total 'gC(v^,/X^) and ^gC(v^,/i+) cross sections 
measured with the recent NOMAD experiment [2]. Due to the high energy of its broad-band neutrino (antineutrino) 
beam with the mean value of 25.9 GeV (17.6 GeV), the normalization was precisely determined from a large sample of 
the deep-inelastic-scattering and inverse-muon-decay events. Using a drift-chamber detector, a total of 14021 (10358 
single track and 3663 double track) neutrino and 2237 antineutrino events surviving the CCQE cuts were recorded. 

In NOMAD, the neutrino events with more than two tracks (the charged lepton plus one proton with momentum 
> 300 MeV) detected would have been removed from the CCQE event sample. On the other hand, in MiniBooNE, all 
events with no pions detected, including two (or more) nucleon knockout events, contribute to the CCQE cross section. 

Additional nucleons may originate from interactions between the struck nucleon and the spectator system in the 
final state, from the nucleon-nucleon correlations in the initial state, and from reaction mechanisms other than nucleon 
knockout, such as those involving meson-exchange currents. 

The effect of multiproton states due to FSI on the CCQE event sample was estimated in Monte Carlo simulations of 
the NOMAD experiment to be very small [2]. Because the nucleon-nucleon correlations produce in nuclei overwhelm- 
ingly proton-neutron pairs, they may have an influence on the track-based measurement of the CCQE cross section 
of neutrinos, but (in the absence of FSI) not on that of antineutrinos. Analyzing the nucleons with momentum larger 



than 300 MeV, we find that in the kinematics of by the NOMAD experiment, scattering off strongly correlated nucleon 
pairs yields ~6% of the inclusive cross sections. This number should be considered as an upper bound of a difference 
between the NOMAD- and MiniBooNE-reported cross sections that the approach of this paper is able to describe. 

Figure 2(b) shows that our calculations of the total inclusive CCQE v^ and v^j cross sections are in good agreement 
with the data obtained by the NOMAD Collaboration [2]. We want to emphasize that their normalization is not scaled 
by additional factors. Although the SF results are higher by ~6% than the NOMAD best fit, this difference is less that 
the ^^% (~11%) systematic uncertainty of neutrino (antineutrino) data. We observe that subtraction of the correlated 
contribution would bring the neutrino result in perfect agreement with the experimental points. However, this may be 
a pure coincidence, as the analogical difference for antineutrinos cannot be explained in the same manner Additionally, 
the results of Ref. [13] suggest that the approach of this paper may, to some extent, overestimate the cross section, 
owing to inappropriate description of the low- 2^ contribution. 

From the total CCQE v^ cross section, the NOMAD Collaboration extracted the axial mass 1.05 ±0.02(stat)± 
0.06(syst) GeV. We estimate that in the SF approach that would correspond to 1.17 GeV, due to stronger quenching 
of the cross section. This value has been obtained without subtracting the correlated strength and, as such, gives the 
lowest Ma required to fit the NOMAD data using our approach. 

Analyzing the shape of the g^g^, distribution of double-track CCQE v^ events for 0.2 < g^g^, < 4 GeV^, the NOMAD 
Collaboration found M^ = 1 .07 ± 0.06(stat) ± 0.07(syst) GeV. Within the quoted uncertainties, this value is marginally 
consistent with our estimate based on the total cross section. 



DISCUSSION 

The axial form factor Fa{Q^) is typically parametrized in the dipole form. Its value at Q^ ^0 is known from the 
neutron beta-decay measurements, whereas the dependence on Q^ is governed by the axial mass M^- 

Because of the sizable contribution of the axial form factor to the (anti)neutrino cross section, the axial mass may 
be determined from the total cross sections or from shape fit to the event distribution with respect to some kinematic 
variable, e.g., Q^. For free nucleons these two methods are equivalent, provided the dipole parametrization holds true. 
For bound nucleons, however, this may no longer be the case, due to inaccuracies of the applied description of nuclear 
effects. Note that different approaches could yield nearly identical (single) differential cross sections and, at the same 
time, different values of the total cross section, or vice versa. 

We made use of this observation, showing for CCQE interaction that, within the SF approach, the axial mass 
extracted by the MiniBooNE Collaboration from the first shape analysis of the g^g^, event distribution [13] is in good 
agreement with the total ^gC(v^,/X^) cross section measured with the NOMAD experiment [2]. 

In Ref. [13], the MiniBooNE Collaboration interpreted the higher value of the axial mass, extracted from the shape 
of the g^gg distribution of CCQE(-like) v^ events, as an effective method of accounting for the nuclear effects neglected 
in the REG model. In fact, the processes which cannot be accounted for within the lA framework have recently been 
shown to contribute to CCQE-like scattering, increasing the total cross section and decreasing the slope of the Q^^ 
distribution of events. Nieves et al. [10] have observed that these effects may effectively be described using a higher 
value of the axial mass, confirming the interpretation of the MiniBooNE Collaboration. 

Nevertheless, the results of Nieves et al. [10] suggest that the neutrino flux in the MiniBooNE experiment is 
underestimated by 9 .2 ± 3 .5 % when Ma = 1 .077 GeV is used in the approach based on the local REG model. Because 
the total cross sections obtained in the REG model are typically higher by ^-^10% than those in the SF approach, the 
20% normalization discrepancy that we have observed in comparison to the MiniBooNE data seems to be in perfect 
agreement with the finding of Nieves et al [10]. Note that the MiniBooNE flux estimate is based on the extrapolation 
of the cross section for ;r* production in p-Be scattering measured on a thin target to a target 35 times thicker 
Such calculations are known to involve severe difficulties [15]. The idea that the MiniBooNE flux might have been 
underestimated is not new [14]. Should it turn out to be correct, this would reconcile the MiniBooNE and NOMAD 
data in a simple way. 

In various interaction channels, the ratio of the events recorded with the MiniBooNE experiment to those predicted 
in the Monte Carlo simulation exceeds unity. It is interesting to note that the size of the ratio is similar: 1 .21 ± 0.24 in 
CCQE scattering [13], 1.23 in CC charged-pion production [16], and 1.58 ±0.05(stat)±0.26(syst) in CC neutral-pion 
production [17]. The quoted figures refer to the simulations with Ma — 1.23 GeV and the parameter K, governing 
the enhancement of the Pauli blocking effect, set to 1.019. We checked that our approach and the model applied by 
MiniBooNE yield the total CCQE cross section differing by < 2.7% for the neutrino energy higher than 350 MeV. 
Therefore, we may conclude that our 20% normalization discrepancy is in good agreement with the MiniBooNE 



simulations and does not seem to be limited to NCE and CCQE interactions. 

The latter conclusion is supported by the measurement of the inclusive CC V^ cross section with the SciBooNE 
experiment [18], using the same neutrino beam as MiniBooNE. The obtained result is larger than its Monte Carlo 
estimate by a factor of 1.12 (k ~ 1.0) or 1.29 (k ~ 1.022), depending on details of the simulation. 

Our calculations underestimate the absolute values of the MiniBooNE-reported cross sections, while being in a good 
agreement with their shape and with the NOMAD data. In the NOMAD analysis, CCQE v events are defined as 
containing at most one proton detected. When the proton's kinetic energy is measured to exceed 47 MeV (momentum 
> 300 MeV), the event is classified as a double-track one. The cross section extracted from the single- and double-track 
events, composing 73.9 and 26.1% of the collected sample, respectively, is shown not to differ. Therefore, only those 
two- and multinucleon final states (2NFS) which involve additional protons of kinetic energy lower than 47 MeV and 
any neutrons contribute to the NOMAD result. Moreover, the contributions of the 2NFS to the single- and double-track 
samples seem to be equal and do not show energy dependence. 

In the MiniBooNE analysis, CCQE events are required exclusively not to involve pions. Hence, the cross sections 
may be increased by a broader class of 2NFS, involving two or more protons, each of kinetic energy > 47 MeV. 
As a consequence, the 20% discrepancy between our calculations and the MiniBooNE data could, in principle, be 
ascribable to more sizable contributions of 2NFS than accounted for in our approach by using the effective axial mass 
of 1.23 GeV Lacking apparent dependence on energy, it is, however, constrained by the NOMAD results. Note that 

(i) the discrepancy in the NCE channel stems from the final-state nucleons with 50 < T < 650 MeV, where T is their 

kinetic energy in total, and appears to the same extent over the whole range of T; 
(ii) because the contributions of 2NFS to the MiniBooNE-reported NCE and CCQE cross sections seem to be equal, 

the kinematics of the knocked-out nucleons in these two cases should not differ significantly; and 
(iii) the NOMAD results constrain the missing strength to involve at least two protons with T > 94 MeV. 

The discrepancy between our calculations and the NCE cross section measured with MiniBooNE is not limited to 
r > 94 MeV, corresponding to Q^ > 0.177 GeV^ and remains constant at the interval 50 < T < 650 MeV, which 
suggests that all the contributing 2NFS channels are open below T ^ 50 MeV. These features do not seem to be 
consistent with the hypothesis that 2NFS contribute to the MiniBooNE and NOMAD data in a different manner, but 
rather point to the flux uncertainty in MiniBooNE being higher than reported. 

SUMMARY 

In this paper, we have applied the spectral function approach to describe nuclear effects in (anti)neutrino scattering 
off carbon nucleus, treating in a consistent manner NCE and CCQE interactions. We have considered a broad energy 
range, from a few hundreds of MeV to 100 GeV. The dipole parametrization of the axial form factor with the cutoff 
mass 1 .23 GeV has been used, as determined from the shape of the Q^^ distribution of CCQE v^ events by the 
MiniBooNE Collaboration in Ref. [13]. This effective method of accounting for two- and multinucleon final-state 
contributions to the cross section seems to be justified in view of recent results of Nieves et al. [10]. 

It has been shown that our approach provides a fairly good description of the NCE v and v differential cross 
sections d<j/dQ^ measured with BNL E734. A good agreement has been found with the total CCQE v^ and v^j cross 
sections from NOMAD. While our calculations provide very accurate description of the shape of the NCE neutrino 
differential cross section da/dQ^ obtained from MiniBooNE, they underestimate its absolute value by 20%. The same 
discrepancy is observed for the flux-unfolded total CCQE cross section reported by the MiniBooNE Collaboration. 

The difference between the total 'gC(V|j,jU^) cross section from NOMAD and MiniBooNE is sometimes attributed 
to a sizable contribution of two- and multinucleon final states. This reasoning is based on the fact that in the latter 
experiment, nucleons knocked out from the nucleus in CCQE neutrino scattering have not been detected. 

We have argued, however, that the NCE data provide an indication of the kinematics of nucleons in CCQE scattering, 
because the MiniBooNE results strongly constrain the allowed differences between nuclear effects in NCE and CCQE 
Vfx interactions. In the kinematic region used to extract the NCE cross section, we find no evidence for the contribution 
of two- and multinucleon final states which would not have contributed also to the CCQE cross sections reported 
by the NOMAD Collaboration. Therefore, the discrepancy between the results from the MiniBooNE and NOMAD 
experiments seems more likely to be ascribable to underestimated flux uncertainty in the MiniBooNE data analysis. 



ACKNOWLEDGMENTS 

The author would like to thank Omar Benhar for illuminating discussions and careful reading of the manuscript. 
Thanks are also due to the Organizers of Nulntl2 for providing partial financial support. This work was supported by 
the INFN under Grant No. MB3 1 . 



REFERENCES 

1. A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev. D 81, 092005 (2010). 

2. V. V. Lyubushkin et al. (NOMAD Collaboration), Eun Phys. J. C 63, 355 (2009). 

3. E. F. Martinez and D. Meloni, Phys. Lett. B 697, 477 (2011). 

4. A. M. Ankowski, O. Benhar, and N. Farina, Phys. Rev. D 82, 013002 (2010). 

5. K. Abe etal, Phys. Rev. Lett. 56, 1107 (1986); 56, 1883(E) (1986); L. A. Ahrens etal, Phys. Rev D 35, 785 (1987). 

6. A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev D 82, 092005 (2010). 

7. A. M. Ankowski, Phys. Rev C 86, 024616 (2012). 

8. O. Benhar, A. Fabrocini, S. Fantoni, and I. Sick, Nucl. Phys. A 579, 493 (1994). 

9. O. Benhar, A. Fabrocini, and S. Fantoni, Nucl. Phys. A 505, 267 (1989). 

10. J. Nieves, I. R. Simo, and M. J. Vicente Vacas, Phys. Lett. B 707, 72 (2012). 

11. O. Benhar, N. Farina, H. Nakamura, M. Sakuda, and R. Seki, Phys. Rev D 72, 053005 (2005). 

12. A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev D 79, 072002 (2009). 

13. A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev Lett. 100, 032301 (2008). 

14. S. Kopp, talk given at the Sixth International Workshop on Neutrino-Nucleus Interactions in the Few-GeV Region (NUINT09), 
Sitges, Spain, 18-22 May 2009; the presentation available at http : //nuint09 . if ae . es/Agenda . html. 

15. M. ApoUonio et al. (HARP Collaboration), Phys. Rev C 80, 065204 (2009). 

16. A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev D 83, 052007 (201 1). 

17. A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev D 83, 052009 (201 1). 

18. Y. Nakajima et al. (SciBooNE Collaboration), Phys. Rev D 83, 012005 (201 1). 



